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We consider a mixture of two-component Fermi gases at low temperature. The density profile of 
this degenerate Fermi gas is calculated under the semiclassical approximation. The results show that 
the fermion-fermion interactions make a large correction to the density profile at low temperature. 
The phase separation of such a mixture is also discussed for both attractive and repulsive interatomic 
interactions, and the numerical calculations demonstrate the exist of a stable temperature region 
T c i < T < T C 2 for the mixture. In addition, we give the critical temperature of the BCS-type 
transition in this system beyond the semiclassical approximation. 
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INTRODUCTION 

Since the experimental realization of Bosc-Einstein 
condensation in dilute gases of rubidium[l-4], 
sodium[5,6], lithium[7], and hydrogen[8], a great 
deal of interest in trapped ultra-cold atoms has con- 
centrated on the topic of trapped degenerate Fermi 
gas. However, it is difficult to achieve a degenerate 
state for fermonic atoms, because the s— wave collisions 
between fermions in a same state are suppressed by the 
Pauli principle, and the p— wave scattering as well as 
the dipole-dipole magnetic interaction are very weak 
at low temperature. The successful demonstration 
of overlapping condensates in different spin states of 
rubidium [9,10] and sodium [11] open a door to study 
the degenerate fermionic gas, because we can cool down 
one component of a mixture by sympathetic cooling. 
Inspired by this observation a number of experiments 
has been conducted on systems of Bose-Fermi mixtures, 
and most recently, using two-component evaporative 
cooling strategy, DeMarco and Jin [12] have succeeded in 
cooling fermonic atom gas down to about 0.5 Tp(300nK , 
depends on the trapping frequencies and the number of 
the trapped atoms). Below this temperature quantum 
degeneracy behaves as a barrier to evaporative cooling 
and as a modification of the classical thermodynamics. 
For the experiment reported by Demarco and Jin[12], 
the atom 40 K were trapped in two magnetic sublevels, 
\F = 9/2, Trip = 9/2} and [9/2,7/2), this mixture of 
atoms states is metastable against mp which changes 
collisions at low temperature, therefore the atom in each 
state is separately conserved. In such a two-component 
mixture of trapped spin-polarized K atoms, interac- 
tions between atoms in different hyperfine states are 
much larger than those among atoms in the same state. 
Indeed, under this approximations a relatively high 
temperature T c for a BCS-type phase transition was 
predicted[13-14]. 



The purpose of this paper is to examine the properties 
including normal and BCS-type phase transition of such 
trapped atoms. The normal state properties of such a 
system were studied under the semiclassical approxima- 
tion in Ref. [15-18]. However, as the interactions between 
40 K atoms in the two different hyperfine states are con- 
siderable strong, it is important to include the effect of 
these interactions in any realistic treatment of the sys- 
tem. The present paper extends the analysis of Ref. [13- 
18] by considering both the discrete nature of the trapped 
atom and the effects of the interactions among them. 

The paper is organized as follows. In Sec. II, we 
analyze the influence of the trap potential and the 
interactions on the density profile of the trapped atomic 
gas. We show that the cloud of the trapped atoms is 
compressed (diluted) for the case of attractive(repulsive) 
interactions. The stability properties of the trapped 
two-component Fermi gases are considered in Sec. Ill, 
and in Sec. IV we investigate the BCS-type transition in 
the system by taking the discrete nature into account. 
We find that the discrete nature is indeed make sense, 
they decrease the BCS-type transition temperature. 
The results for the BCS-type transition are beyond the 
semiclassical approximation. Finally, we summarize our 
results in Sec. V. 



DENSITY PROFILE OF A TRAPPED 
INTERACTING TWO-COMPONENT FERMI 
GAS 



We consider a dilute gas which consists of interacting 
fermionic two-level atoms trapped in an external poten- 
tial Vo(r). As the gas is dilute, the interactions mainly 
happen through two-body collisions. Furthermore, be- 
cause the s— wave scattering length between fermions in 
a state is suppressed, and the p-wave scattering is greatly 



2 



reduced due to the presence of the centrifugal barrier, 
we may neglect the interactions between fermions in the 
same hyperfine state and only consider the s-wave scat- 
tering between the fermions in different hyperfine states. 
Under this consideration, the system is then described by 
the following Hamiltonian 

H = Hi + H 2 + Hi n t , 

Hi = J dr4{^- ii + V {r))^{r),{i = 1,2) 

H int = v J drdr ip\(r)tpl(r' )S(r - r' )ip 2 (r' )ipi(r){l) 

where m is the mass of one fermion, and the interatomic 
potential has been approximated by a constant poten- 
tial Vo5(r — r ). ipi(r) stands for the annihilation opera- 
tor of a fermion at position r in the hyperfine state \i), 
and it obeys the usual fermionic anticommutation rules. 
The trapping potential is for simplicity taken to be an 
isotropic harmonic oscillator potential Vo(r) = imu> 2 r 2 , 
and the trapping frequency is the same for each hyper- 
fine state. In addition to what we stated above, we have 
assumed that the number of particles N in each state is 
the same such that we only have one chemical potential 
/i. As the critical temperature for a BCS-type transition 
is maximum when the number of particles in the two 
hyperfine state is equal[13,14], so this configuration has 
most experimental relevance. Indeed, in the experiment 
reported in Ref.[12], the number of the atoms in state 
mp = 9/2 and mp = 7/2 are approximately equal. The 
noninteracting case is achieved by setting v n — 0; this 
limit has been discussed in ref.[15,16,18] for one compo- 
nent Fermi gas within the scmiclassical approximation. 
In this section, we are interested in the effect of the inter- 
action and the trapping potential on the density profile 
under the semiclassical approximation. We can, there- 
fore, ignore any pairing correlations leading to BCS-type 
transition, and use the mean field Hamiltonian 

H m i = - ^+ -muj 2 r 2 ) +v (iPl(r)^ 2 (r)), 

H m2 = (^-»+^mu; 2 r 2 )+v (4(r)Mr)), (2) 

this equation comes from eq.(l) straightforwardly. Here, 
H m i describes the effective Hamiltonian for component 
i. (ipl (r)ipi(r)) is the standard Hartree-Fock result for 
a hard sphere interaction model. In order to get some 
analytical results, we study the density profile here within 
the semiclassical approximation, it is given from eq.(2) 
that 

n i (r) = ^(^)l/3- | /|(^)(i = l,2), (3) 

where (3 — l/fc^T, ks is Boltzmann's constant. z\ = 
exp(l3(ii-v n 2 {r)-V a {r))), and z 2 = exp((3(^-v ni(r)- 
Vb(r))) are called local fugacity[19] for component 1 and 



2, respectively. the chemical potential for the atom, is 
determined through Ni — J drrii(r). Before we present 
the density profile of the two-component Fermi gases, we 
consider a range of parameters for relevant experiment 
reported in Ref. [12]. The potentials for the centre-of- 
mass motion of a single atom in the hyperfine state can 
be approximated as a cylindrically symmetric harmonic 
potential with an axial frequency of lo z — 2ir x 19.5-ffz 
and a variable radial frequency, which can be varied from 
oj r = 2tt x MHz to 2n x 370Hz. With the temperature 
being cooled down, the quantum statistical properties of 
the trapped gases become more evident, and at temper- 
ature T ~ 0.5T F ~ 300nK, the effect of Fermi-Dirac 
statistical are observed in the momentum distribution 
of the gas. With these parameters, it is evident that 
fiw/ksT « 1, i.e., the semiclassical approximation is a 
good approach to the realistic case. (This does not in- 
dicate that the semiclassical approximation holds well at 
the temperature where BCS-type transition occurs) As 
the gas is dilute, we can expand rii(r) up to first order of 
the coupling constant vo, one gives 

n a (r) = ^(-^)i / 3- | (/|(^°)--o(^ ) 2 / i (^)K), 

(4) 

where 

^(r) = |(^)^- | /|(^°), 
2 nn 2 

and a,b = 1,2, z° = z° = exp{{3{pL - \muj 2 r 2 )). The 
result shows that for i; < 0, i.e., attractive interaction, 
the cloud of particles is compressed as compared to the 
noninteracting result. Because a high density of particles 
increases the critical temperature for a BCS-type tran- 
sition, this effect favors the formation of the superfluid 
state [13]. 



STABILITY PROPERTIES OF A TRAPPED 
TWO-COMPONENT FERMI GAS 

Since the experimental realization of the two- 
component condensate. Most of the theoretical works 
concerning multicomponent condensates [20-26] has 
been devoted to systems of two Bose condensates and 
Bose- Fermi mixture[27-30]. However, other systems are 
of fundamental interest, one of these is two-component 
trapped Fermi gas. In fact, by using the sympathetic 
cooling, the fermionic atoms in different hyperfine states 
has been cooled down to 300nK, and the quantum 
statistical effect in this system has been reported[12]. 

For the two-component fermion system, the thcrmody- 
namical properties are trivial if there are not interaction 
between them. But in this case the sympathetic cooling 
scheme can not make effect and the degenerate fermions 
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in a trapped potential can not been achieved. The ther- 
modynamical properties may be changed when the inter- 
actions within and between the two hyperfine levels turn 
on. Then a new phenomenon, the phase separation, may 
occur in this system. For a homogeneous fermion mix- 
ture system, the Helmholts free energy can be written 
as[31] 



2V 

(3F = -- T f^(z 2 ) + 2a 2 N 2 p 2 X 2 2 

A 2 



2V 



^3-/ f ( Zl ) + 2a lPl N 1 \( + a 12 (X( + \ 2 )N 2 Ni 



where index 2 refers to the fermionic component in level 
9/2, 7/2), whereas index 1 stands for them in 9/2, 9/2), 
Ni is the number of particles in component i, A denotes 
the thermal wave length of the atoms, /„(z) represents 
the Fermi integral, a\ and denote the coupling con- 
stants. From eq.(5) we obtain the chemical potential for 
each component straightforwardly, 

0m = l3^ + Aa lPl \l + a 12 {\l + \ 2 2 )N 2 /V, 

/3 M2 = /3 Ai o + 4a 2 p 2 A2 + a 12 (A2 + A2)iV 1 /y, (6) 

where $ are the chemical potentials of ideal gas. There 
are three terms in each chemical potential(6). The second 
term comes from the interaction within the component, 
while the third term is from the interaction between the 
fcrmions in different levels. As known, an homogenous 
binary mixture is stable only when the symmetric matrix 
ft given by 



A* = 



dpi 

dpi 
dpi 



dpi 
dp2 
dpi 
dpi 



(7) 



is non-negatively definite. In other words, all eigenvalues 
of matrix p given in eq.(7) are non-negative. Mathe- 
matically, for homogenous two-component fcrmions the 
stability conditions are 



dui dn 2 
dpi dp 2 



and 



det 



dpi dpi 

dpi dpi 

dpi dpi 

dpi dpi 



> 0. 



(8) 



(9) 



For ideal gas, we have p\ = yjr/a (^i), and p 2 = T7/3 (z 2 ), 
this leads to 



A 



dp 2 fi&Y dp 2 /i(zi)' 
It follows from eqs (8) and (9) that 
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A(*i) 



>o, 



(10) 
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Temperature T 



FIG. 1: Z(T, 012) vs. temperature T. The coupling constant 
CI12 is different for Fig.l-a and Fig.l-b. a:ai2 = 0.5, b:ai2 = 
0.2. The parameters for Fig.lc are the same as those in Fig.l- 
b, but for the inhomogeneous case. 



4a 2 A^ 
and 

Z(T,a 12 ) = (4 ftl A2 + 



A? 



2\2 



>o, 



)(4a 2 A2 + 



(12) 



Ai 



ai 2 (\t + y 2 y>o. 



(13) 



Now, we discuss the stability of this system with repulsive 
interactions. The case with attractive interactions will be 
discussed in the next section. It is obvious that the sta- 
bility condition (11) and (12) hold always for a\ > 0, 
a 2 > 0. We would like to point out that the stabil- 
ity conditions (11-13) do not involve the densities of the 
both components. At hrst glance, this seems to be confu- 
sion. In fact, there is no contradiction. One can demon- 
strate that at low density the Hclmholtz free energy of the 
gas reduce to a quadratic form in N\ and N 2 . To have 
a minimum, this form should be positive definite, i.e., 

de^ll dNidNi H — 0- Therefore, the corresponding stabil- 
ity criterion involves only density-independent constants. 
This criterion is similar to the stability conditions for 
two-component Bose-Einstein condensate in a trapped 
ultra-cold gas[23,24,33-37]. When T -> 00, A 4 -> 0, hence 
Z{T,a\ 2 ) ~ ^jj- Thus at high temperature, the homo- 
geneous binary gas mixture is always stable and no phase 
separation occur. The quantity Z{T,a\ 2 ) as a function 
of the temperature is illustrated in Fig.l. We see from 
Fig.l that the system is always stable when T — > and 
T — > 00, and the system is unstable for T c \ < T < T c2 . 
In particular, T c \ and T c2 depend on ai 2 , the scattering 
length for fcrmions in different states. As a\ 2 decreases, 
T c \ tends to T c2 (going from fig.l-a to fig.l-b). 

/ We would like to note that, the stability condition 
depends only on |ai 2 |, so the fermions with ai 2 and — a± 2 
have the same stability condition. And the stability dis- 
cussed above is only a result of interactions among dif- 
ferent components. 
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Until now, we considered only a homogeneous Fermi 
gas mixture at finite temperature. In practise, how- 
ever, experiments with ultracold atoms are performed 
by trapping and cooling atoms in an external potential 
that can be generally modeled by an isotropic harmonic 
oscillator V(r) = ^w 2 r 2 , where ui is the trapping fre- 
quency. An exact criterion for the stability of an inhomo- 
geneous Fermi-Fermi mixture should involve calculating 
the Helmholtz free energy. Fortunately, in the system 
considered here it is a good approach to take use of the 
semiclassical approximation, which treats the atoms as a 
local homogeneous system. This approximation requires 
that the level spacing Ttw of the trapping potential is 
much smaller than the Fermi energy. Of course, the sem- 
icalssical approximation always breaks down at the edge 
of the gas cloud where the density vanishes and the effec- 
tive Fermi energy becomes zero. In this approximation, 
the stability conditions can still be calculated by means 
discussed above, with the understanding that the effec- 
tive chemical potentials arc spatially dependent through 



1 2 2 

Mi = Mi - 2 muJ r ' ^ 2 



M2 



1 



2 2 

muo r . 



In this sense, the stability condition is the same as given 
in eq.(13) but replacing Zi(i — 1, 2) by 



zi = zie 



and z 2 = z 2 e 2 ' 



As shown in Fig.l-c, the region of temperature in which 
the system is unstable becomes narrow as compared to 
the homogeneous case(r = 0). 



BCS-TYPE TRANSITION IN TRAPPED 
TWO-COMPONENT FERMI GAS 

The achievement of atomic Bosc-Einstein condensation 
has induced an experimental growth of interest in the 
properties of ultracold dilute quantum gases. Of partic- 
ular interest now is the physics of trapping and cooling 
of fermionic atoms. Indeed the prospect of superfluid- 
ity with dilute atomic vapors has already been studied 
within the semiclassical approximation by several groups 
[13,14,38,39]. Because the semiclassical approximation 
is not of fundamental quantum physics, we will extend 
the analysis under the semiclassical approximation by in- 
cluding the discrete nature of the trapped atom in this 
section. 

Let us consider two species of fermions in a trap, 
which interact with each other by two-body collisions ( 
s— wave scattering). The Hamiltonian describing such a 
system is given by eq.(l). The two species of fermions 
correspond to the trapped atoms in two hyperfine lev- 
els |1) = |9/2,9/2) and |2) = |9/2,7/2). Expanding 
V£(r)(a = 1,2) by 



(14) 



where creates one particle in state \an) = \a) (gi \n), 
which satisfies 

h 2 1 

(- — V 2 + -mLj 2 r 2 )\an)=e na \an), (15) 
the Hamiltonian in eq.(l) becomes 



where V; 



ijmn 



v J drdr 5(r-r )4>*(r)(j)*(r )(j> m {r )<j> n {r) 
and a, —a denote the two hyperfine levels(say, for ex- 
ample, if a = |9/2, 9/2) , then -a = 9/2, 7/2)). The 
next step in a mean-field treatment of the Hamiltonian 
in eq.(16) is to develop the operator product VanV'-cm 
around their mean values by substituting 

To first order in the fluctuations, we are left with the 
effective mean-field Hamiltonian 

H — ^ ^ £anQa n Q j an~\~~ } ^ ^m a lm fl -affl^ } ] A m a am (2_ am 
an m m 

(17) 

Here A m = (A TO )* = ^2 n V mrnrm {aana—an) is the equi- 
librium value of the BCS order parameter. As the ef- 
fective mean-field Hamiltonian in terms of the operators 
ajj n and a an is non-diagonal, one can not directly cal- 



culate the expection value (aj 



>). This is, as usual, 



resolved by first applying a Bogoliubov transformation 



bam u m&am ~F ^m a — am^ 



bam — U m a lm "F V m a 



m^ — ami 



(18) 



to diagonalize the Hamiltonian in eq.(17). After per- 
forming this unitary transformation, we require that the 
Hamiltonian in terms of the new quasiparticlc operators 
b am and b\ m has only diagonal elements, and further- 
more we assume that these operators obey the usual an- 
ticommutation relations still. This determines the val- 
ues of the yet unknown u m and v m . The latter con- 
straint requires that the constant u m and v m must sat- 



isfies the relations \u r , 



= 1 and the require- 



ment of diagonality of the Hamiltonian after Bogoluibov 
transformation lead to \v m \ 2 = 1 _ em / Em with 

E m = yj A 2 ^ + e 2 ^. E m are eigenvalues of the Bogoluibov 
quasiparticles. 

Using Bogoluibov transformation(18), the equilibrium 
value of the BCS order parameter is calculated easily, it 
is given that 



A n — ^ ^ Vrnmnn tCLTih En 



(19) 



As usual, the order parameter does not depend on its 
index m. So we arrive at the gap equation 



5> 



i 



ztanh 



V c r, 



A 2 



2k B T 



= 1. 



(20) 
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Setting A = 0, one finds the critical temperature T c as 
a function of the trapping frequency a; [18]. The density 
N(0) of atoms near the Fermi surface and the coupling 
constant Vq = V w 



v mmnn : 



v N(Q)ln(1.13/3 c nu) + v Q Bk B T c 
1 



exp[ 



k B T a 



-v Q Bk B T c - v Bk B T c ln(l 



1 

exp[ 



htu ,. 



knT r 



VQBkBT c ln2, 



(21) 



where B = 



and f3 c — 



The first term in the 



2(huj) 2 > Jc k B T c 

right hand side of eq.(21) is just from the usual BCS 
theory, in other words, if the semiclassical approximation 
is a good approach to the theory or the discrete nature 
of the trap levels can be neglected, T c = |^ e _ ^o«(o) _ The 
rest terms in the right hand side of eq.(21) are corrections 
of the discrete trap levels to the usual BCS theory. 

This superfluid phase transition, which is similar to 
the BCS transition in a superconductor, might occur at 
very low temperature. At such low temperature, whether 
the semiclassical approximation holds depend on both 
the temperature and the trapping frequency. Hence 
we investigate here the BCS-type transition from the 
other aspect, beyond the semiclassical approximation. 
The results show that if the semiclassical approximation 
holds, i.e., huj/kBT c << 1, the transition temperature is 
just the BCS one. Otherwise the effects of the discrete 
trap levels provide a negative correction to the transition 
temperature. 



CONCLUSION 

In summary, we considered a Fermi gas occupying 
two hyperfine states trapped in a magnetic field. Atoms 
in different hyperfine levels can interact via s— wave 
scattering. Under the semiclassical approximation, 
we calculated and discussed the density profile for 
the trapped fermions. The purterbative results up to 
the first order of the coupling constant show that the 
density of the atom is compressed(diluted) as compared 
to the noninteracting case due to attractive (repul- 
sive)interaction. We also investigate the mechanical 
and statistical stability of the two-component gas with 
interaction between the atoms in different hyperfine 
levels, and find that these interactions strongly affect 
the stability of the system at finite temperature. The 
regime T c % < T < T c i in which the system is unstable 
depend on the strength of the interaction and the spatial 
atomic position. Furthermore, we consider the BCS-type 
transition beyond the semiclassical approximation, 
within which the most current literature study the 
superfluid state of trapped fermonic atoms. The results 



showed the transition temperature is a function of 
the trapped frequency, the coupling constant and, as 
usual, the density of the atoms near the Fermi surface. 
Neglecting the discrete nature of the trap levels, the 
transition temperature return to the BCS's one, and the 
correction to T c is negative due to the discrete nature of 
the trapped atom. 
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